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In this paper, we study a number of stochastic models including station- 
ary, nonstationary, and linear processes for the purpose of simulating 
earthquake- or explosion-induced ground motions. The important sta- 
tistical characteristics of each model and their effects on structural systems 
are investigated in soyne detail. We obtain expressions for the mean-square 
response of simple linear mass-spring oscillators to each model. We discuss 
numerical procedures for time series simulations of these models. The 
objectives of this paper are (i) to examine and compare the statistical 
properties and effects on structures of all possible stochastic processes 
applicable to model ground motions and (ii) to offer engineers a basis for 
forming their own judgments. 

I. INTRODUCTION 

A primary concern in problems dealing with earthquake or blast 
response is the proper definition of random force environments. Typical 
problems include the design of earthquake-resistant frames for electronic 
or mechanical facilities in a building, and the estimation of structural 
damage resulting from nuclear detonation in a given area. Normally 
it is necessary to create ground-motion data artificially from information 
derived from limited recordings. Therefore, in ground-motion analysis, 
as in many other fields of engineering physics, time series modeling or 
simulation problems emerge. This paper carries forward the concern 
with such problems. 

In a previous paper the statistical characteristics of a collection of 
earthquake ground-motion data were analyzed. 1 It was shown that a 
stationary random process of finite duration could be used to model 
the high-intensity phase of a ground-motion accelerogram. It was sug- 
gested that a narrowband stationary process be used. However, the 
narrowband process, like many other stationary models, fails to produce 
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the initial buildup and the decaying terminations which are apparent in 
many real ground-motion accelerograms. To simulate these portions, 
nonstationary models must be used. In this paper we investigate a 
number of physically realizable stationary, nonstationary, and linear 
stochastic models, when applicable to the simulation of ground motions. 
We derive the important statistics of these models and discuss numerical 
simulation procedures. The covariance, autocorrelation function, and 
power spectral density of the response process of a class of linear time- 
invariant systems excited by the input process represented by each 
model are examined. It is hoped that by comparing these models for 
basic definition, for simulation procedure, and for effects on induced 
structural responses, a practicing engineer will be able to select one that 
will be suitable for the analysis of his specific problem. It is the intent 
of this paper not only to sum up the current state of the art/science but 
also to give to researchers in the field of ground motion study such hints 
and directions for the development of advanced stochastic models as 
mathematical sophistication permits. 

II. STOCHASTIC MODELS 

2.1 Stationary Models 

A typical ground-motion accelerogram from an earthquake or nuclear 
detonation consists of three phases: an initial rapid-buildup phase, a 
high-intensity primary phase, and a gradually decaying tail — all ob- 
viously nonstationary phenomena. However, many researchers still 
prefer to model earthquakes by stationary processes because the low- 
amplitude starting and ending portions of an accelerogram do not sig- 
nificantly affect the structural response as compared with the response 
induced by the primary phase. Therefore it appears legitimate to model 
the primary phase by a stationary process of finite duration. Further- 
more, the stationariness assumption greatly simplifies the response 
evaluations and numerical simulation procedures. This feature is partic- 
ularly important from the practical viewpoint. Three stationary models, 
designated as x t (t), x 2 (t), and x 3 (t), along with their autocorrelation 
function and power spectral density are defined in Table I. 

The first model, white noise Xi(t), defined as a stationary random 
signal having gaussian probability amplitude distribution and a con- 
stant spectral density for all frequencies, is the simplest one of all. 
Numerically, it can be simulated by generating a sequence of gaussian 
independent samples of gaussian random numbers g n , spacing them at 
small time interval At, and assuming linear variation between amplitudes 
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over each At. Bycroft 2 has studied this overly simplified model on an 
analog computer. This model fails to provide any frequency descriptions 
of the motions that are so important to the structural response analysis; 
its use results solely from its mathematical simplicity. 

For model x 2 (t), the filtered white noise, symbols h (t) and H n (u) 
in Table I represent the transfer function in the time and frequency 
domains, respectively. The process x 2 (t) is a gaussian, co variance sta- 
tionary, narrowband process. The use of x 2 (t) in modeling earthquake 
ground motions is based upon the resemblance of the autocorrelation 
function, the power spectral density, and the response spectra of strong 
motion earthquakes to those of the narrowband process. 1 ' 3,4 The 
numerical simulation of x 2 (t) based on prescribed power spectral density 
is a routine exercise using the basic relationship between the power 
spectral densities S x i (w) and S x2 (u) of the input and output of the linear 
filter: 

&,(«) = \H (u)\ 2 S xl (o>). (1) 

This equation suggests that a signal x 2 (t) can be created from the random 
process Xi(t) whose spectral density S x i(ui) = 1 by passing x,(t) through 
a filter whose transfer function H (co) satisfies | H (u) | 2 = S x2 (ui). 
A detailed approach to simulating stationary processes with a rational 
spectral density function, which is represented by a quotient of two 
polynomials in w, has been described by Franklin. 8 Methods of estimat- 
ing ground-motion spectral densities are given by Liu and Jhaveri. 1,4 
A sample function of process x 2 (t), generated digitally by passing a white 
noise through a linear filter with a natural frequency of 10 rad/s and 
a damping ratio of 5 percent, is shown in Fig. 1 which demonstrates a 
typical appearance of a narrowband process. The justification for using 
this process to model ground motions is that both its power spectral 
density function and its response spectrum resemble and can be made 
to match those of real-world records. It reflects the predominant effects 
of the site on motion and, by using a modal analysis, x 2 {t) can model 
ground motions which exhibit multiple peaks in their frequency spec- 
trograms. 4 

For x 3 (t) as defined in Table I the random impulses \a n ] form a 
stationary discrete parameter process, and {t n } is a stationary point 
process independent of |a„}. 6 ' 7 As in the expressions for autocorrelation 
function and power spectral density, /3 is the average number of impulses 
per unit time interval and p n (t) represents the probability density 
function for n consecutive intervals of {t n } within a time duration t; 
p{n) and P„(iw) are defined as 
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p(n) = E[a m a m+ „] (any m) 
P n (iw) = [ p„(0 exp {-tat) dt 



(2) 
(3) 



in which E[ ] denotes ensemble averages. 

Notice that the model x 3 (t) allows correlation p(n) to exist among 
the random amplitudes a n which are assumed to be independent for 
many other stochastic models. The simulation of x 3 (t) can be achieved 
by the spectral approach as described by Franklin or by the correlation 
approach based on the matrix factorization procedure as proposed by 
Moore and Anderson. 8 However, when p n (t) is poisson with mean arrival 
rate and Markov correlation p(n) - p |n| and \p\ ^ 1, a sample func- 
tion of x 3 {t) can be created by simultaneous simulations of independent 
samples of \a n ] and {t n \. In this case the probability density function 
of the waiting time r k = fc+i — t* is an exponential distribution, that is, 
p( rt ) = ^ exp (— /3r fc ) and therefore a sequence r k can be generated 
from a sequence of uniformly distributed random numbers w k : 



r k = -- In (1 - w k ), 



S -Wk < 1- 



(4) 



The correlation p |n| can be introduced into \a n \ by using an autoregres- 
sive transformation of a sequence of gaussian random numbers g n with 
zero mean and unit variance: 

«, = d- tM (5) 

a' n = pa' n - x + g„ . 




12 16 20 

TIME IN SECONDS 



Fig. 1 — Sample function of process x 2 (t). 
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A digitally simulated sample of x 3 (t) with p = 0.5, /3 = 1.0, and 
poisson distribution for p n (r) is shown in Fig. 2(a); its autocorrelation 
function and power spectral density are shown, respectively, in Figs. 
2(b) and (c). The time response of a linear filter having a natural fre- 
quency co = 10.0 rad/s and a damping ratio f = 5 percent to this input 
sample member of x 3 (t) is shown in Fig. 3(a). The corresponding auto- 
correlation function and power spectral density of the sample response 
are shown in Figs. 3(b) and (c), respectively. Notice that Fig. 3(b) 
exhibits a damped oscillatory motion with a frequency o> = 10.0 rad/s 
which is the natural frequency of the filter. At this frequency there is a 
peak in the power spectral density as shown in Fig. 3(c) which is similar 
to that of strong-motion earthquakes. Notice also that in comparison 
with x 3 (t), which shows abrupt peaks and dips in the waveform, the 
filtered impulse process more closely resembles a ground-motion ac- 
celerogram. Based on these results, it appears that a filtered impulse 
process may be more appropriate than the impulse process x a (t) itself 
in modeling earthquake motions. 

When ground motions are represented by stationary processes, it is 
extremely important to properly determine the duration and intensity 
of the processes because the induced response of structures depends 
heavily on these two parameters. To those who are reluctant to neglect 
the nonstationary effect resulting from the starting and tail portions of 
earthquake accelerograms, the stationary models I through III ob- 
viously are not satisfactory. One might therefore consider using the 
following nonstationary models. 

2.2 Nonstationary Models 

Five useful nonstationary models, x 4 (t) through x 8 (t), are defined in 
Table II. 

The model x 4 (t) is a frequency-modulated nonstationary function. 
For its definition in Table II, &,• , a,- , and w f - are given sets of real positive 
numbers, and iv t are independent random variables from uniform dis- 
tribution over and 2ir. The use of x 4 (t) is based on the result that the 
skewed bell shape covariance function of x 4 (t) is similar to and can be 
made to match that of an earthquake. Numerical simulations of x 4 (t) 
present no special difficulty as its sample members can be created di- 
rectly for given sets of 6, , a,- , and co,- according to its definition. It may 
be expected that use of a large number of terms as given by x 4 (t) will 
produce member functions that look much more like real earthquake 
accelerograms. 9 

For the model x 5 (t), <f>(t) is a deterministic or envelope function and 
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Fig. 2(a) — Digital realization of sample x 3 (t), the impulse process. 

Fig. 2(b)— Autocorrelation function of sample x 3 (l). 

Fig. 2(c)— Power spectral density function of sample x 3 (t). 
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Fig. 3(a) — Digital realization of sample response of Xi(t). 

Fig. 3(b)— Autocorrelation function of sample response of x»(t). 

Fig. 3(c)— Power spectral density function of sample response of x 3 (t). 
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f(f) is a stationary random process with given autocorrelation function 
and power spectral density. The function f(t) may be any one of proc- 
esses Xi(t), x 2 (t), and x 3 (t). This model in its various forms has been 
studied by Peterson and Pullen 10 ' 11 and by MacNeal, and others. 1 " 
The process x 5 (t) has been suggested for model earthquake motions by 
Shinozuka and Sato 13 and by Jennings, and others. 14 The process x a (t) 
is the response function of x s (t). In Fig. 4 is shown a sample member 
of x 5 (t) with tf>(t) = sin (vt/30) and /(£) as the sample member of x 2 (t) 
in Fig. 1. Its response function to a linear filter with a natural frequency 
of 10 rad/s and a damping ratio of 5 percent is shown in Fig. 5 which 
represents a sample of the process x a (t). Figure 4 clearly illustrates 
that the waveform is enveloped by a half-sine wave which produces a 
strong phase between the times of 17 and 23 seconds. The high-ampli- 
tude tail observed in Fig. 5 is obviously undesirable in modeling ground 
motion. However, notice that this figure shows only the response history 
cutoff at the end of excitation. If free vibration is allowed after the 
termination of the input, the highly amplified portion will gradually 
decay, and the resulting waveform will then compare more favorably 
with real ground-motion records as shown. 

A nice feature of process x 5 (t) is its separable property which greatly 
simplifies the mathematics required to evaluate the response statistics. 
Furthermore, the envelope function <f>(t) can be chosen so that the pat- 
tern of the rise and fall of the simulated waveforms is similar to that of a 
real earthquake motion. However, because the choice of #(£) is arbitrary, 
when using x 5 (t) the associated nonstationary effects resulting from the 
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Fig. 4 — Sample function of process xt(t). 
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Fig. 5 — Sample function of process x 6 (t). 

total time duration, the time proportion of the three distinct phases of 
the ground motion, and the rate of buildup and decay must be carefully 
examined. 

The shot noise x,(t) and the filtered shot noise x B (t) are the nonsta- 
tionary counterparts of the stationary white noise Xi(t) and filtered 
white noise x 2 (t), respectively. The nonstationariness is introduced into 
these models by two sources: the time-dependent mean arrival rate 
fl(t) and the time-dependent amplitude joint probability density func- 
tion p(cti , U ; a,- , t t ). However, it is assumed that the impulse ampli- 
tudes are completely uncorrelated, that is, the joint probability density 
function is separable: 

p(a, , h ; a 2 , t 2 ) = p(fl! , ti)v( a 2 , k)- 

Further, the probability of n impulses in a small interval dt is negligible 
for n ^ 2. Based upon these conditions, it can be shown that the covari- 
ance functions of x 7 (t) and x 8 (t) are also time-dependent: 



Cov, 8 (/ 1 ,t*) = j h (t x - T)h (t 2 ~ t)E[ci 2 (t)]0(t) dr; 
Cov l7 (*, ,t 2 ) = j 8(1, - r) S(U - T)E[a\r)MT) dr 



(6) 



= 1(1,) 8(1, - t 2 ); 



(7) 



where r denotes the appropriate time domain, and I(t) is the strength 
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function of the shot noise defined as 

E[x 7 (t)x 7 (t + r)] = I(t) «(t). 

Notice that, when I(t) is a constant, both x 7 (t) and x 8 (t) become sta- 
tionary. 

The process x 7 (t) can be constructed by generating independent gaus- 
sian random variables a n with zero mean and variance equal to I(t n ) At, 
and by linearly connecting them over At along the time axis. The 
process x 8 (t) can be obtained in a similar way simply by shaping each 
impulse with the transfer function h (t) according to its definition in 
Table II. Amin and Ang 15 have used process x 8 (t) with a second-order 
filter to model earthquakes. Both x 7 (t) and x 8 (t) are justified on the basis 
of the similarity between real and nonstationary waveforms and on 
the matching (i) of their time-varying covariance functions and (ii) 
of the induced response spectra with those of the real ground-motion 
data. 

2.3 Linear Models 

Special cases of the mixed autoregressive-moving average process 
x t of order (m, n) as defined below also can be used to model ground 
motions: 

x. - E *i*i-i + gt + E Bi9i-4 , * - 0, ±1, ±2, • • • (8) 

1=1 « = 1 

where #,• and 0, are characterization parameters, and process \g t } is a 
white noise. An wi-order autoregressive (ar) process is given by the first 
sum of equation (8), that is, 

x t = E <t>&t-i + 9t (9) 

•=i 

and an n-order moving average (ma) process by the second sum 

n 

si = g t - E QiQt-i ■ (10) 

The ar processes given in equation (9) are particularly useful in time 
series simulations because they are very flexible and can be used to 
model a wide range of real-world random data. For example, the auto- 
correlation function of a second-order ar process, involving only two 
parameters, can produce a wide variety of autocorrelation functions. 
Therefore by matching the ground-motion autocorrelation function of 
the damped oscillatory type, one can estimate parameters </>,- and fit 
the observed motion to an appropriate linear stochastic model. 
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III. RESPONSE ANALYSIS 



Because a ground-motion model when chosen will be used to specify 
vibration environments for structural testing and design, it is important 
to examine its effects on the time and frequency response of some rep- 
resentative systems. In what follows we shall compute and compare the 
mean-square response of a class of second-order linear systems to stochas- 
tic inputs represented by models as described in the previous section. 
The linear systems are characterized by two constant parameters, the 
damping coefficient f and a natural frequency w,, , and have transfer 
functions 

h(t) = e, V (-M ) s . np(f ^ 0; ^^ 

= t < 0; 

where p = w B (l — f 2 )*. 

Let the response function to x t (i) be y^t); i = 1, 2, • ■ • , 8, for systems 
defined by equation (11). The mean-square response is given by 

Xtotitfi = f / Ht - r)h(t - r') gi (r, r') dr dr' , 

i= 1,2, ••• ,8. (12) 
The function ^(t, t') for each model is 
9l = 2irS 5(t - r'), 

<7 2 = ft* 2 (r - r') = 2tt ( h (t)h (t + r - r') dt, 

•'o 

g a - B,a(T - r') = 0^(0) 5(t - t') + E P(n)p»(| r - r' |)J , 

ff4 = i Z) rr'a' exp (— a,(r + t')) cos 01,(7' - t), 

1 

(7 5 = *(t)0(t')J?/(t - r'), 

g«= [ [ K(t - djhoir' - d 2 )<t>(6Mo 2 )R f (e 2 - 0.) tf0, </0 2 , 

•'r -'r 
7 = £[a 2 (r)]/3(T) 5(r - r') = 7(r) 5(t - r'), 

g a = f h n (r - 0)E[a(d)W) dd f h (r' - e)E[a(O)]0(6) dO 

+ f h n (r - d)h (r' - 9)E[a 2 (9W(d) dd. (13) 

•'r 
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The explicit expression for the mean-square response to stationary 
models can be obtained by integrations using equations (11) through 
(13): 

E[y 2 M = zh + rS ° 



E[yl(t)] =* 



2f«» 2pu n 
•(f cos 2pt - (1 - f 2 )' sin 2pt - 1/f) exp (-2&J), (14) 

irS l2 (u n ) 



[ 



2r«„ 
1 -- 6XP ( ~ 2(J " r0 | p 2 + 2(o; n f) 2 sin 2 p/ + a; n pfsin2p/)J- 



V 



(15) 



It is assumed in equation (15) that S x2 (u), the power spectral density 
of x 2 (t), is given, and the main contribution of x 2 (t) to the response 
comes from the region around co = co„ . ia And: 



w0[c 



&[0a(0] = /3Ci + /3 exp (-2fw n f)| c 2 cos 2pt + c 3 sin 2pt 
+ £§ (« - a 2 ) sin 2 pi + & cos 2 pt - £&* 



V 

20 2 p 
py 



7 2p fw n _ 

exp (— (fro n + o)i)[(a + fw„) sin pt -\- p cos p<] (16) 



in which a = (1 - p)/3, M = a 2 + (1 - 2f K , 7 = (a 2 + "D* - 
4a 2 f 2 a> 2 n , and 



7 



2f«* co n / T 4fw" 



w, |_2p 2 V 7 2^ T J 

1 [> P r / „ rt 1 A , a/3 PM \l 

- — «. + of -si H — ) • 

p L 7 iw„ \^ 7 /_ 



c, = - 



Notice that, as time increases, mean-square response to each of three 
stationary models approaches its steady-state level, that is, 

ttS 



\im E[y%t)] = 



2fco* ' 
irS x2 (u n ) 



2f«" ' 



i = 1 



i = 2 



% = 3. 



(17) 
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Let the rate of convergence of the mean-square response be defined as 

\im E[y%t)] -E[y*(t)] 



*■ = 



\im E[y-(t)] 



(18) 



and N = co n t/2ir be the number of cycles of the motion required to reach 
€ { . A lower bound to the estimate of N can be easily established for 



processes x,(t) and x 2 (t) as the following: 

f i + f + rd - ff \ 
L (i - r 2 )< J ' 



47Tf 



o < r < i. 



(19) 



From this it is noted that N is independent of the natural frequency 
co„ of the system but is heavily dependent on the damping ratio f . Figure 
6 shows that systems with high damping, when excited by x,(0 and 
x 2 (t), will approach steady-state level faster than those with low damp- 
ing. A system with 10 percent damping will reach 90 percent of its steady- 
state response in two cycles of motion. 

Because the mean-square response to x 3 (t) involves more parameters 
than that to x t (t) or x 2 (t), the response convergence rate for x 3 (t) is 
more difficult to evaluate. However, f remains as the dominant factor, 
not or p. The mean-square response of x 3 (t), as of x^t) and x 2 (t), ap- 
proaches its steady-state value rapidly when f is high, slowly when f is 
low. 

Explicit expression for the mean-square response to nonstationary 
models can also be obtained similarly from equation 13 although the 
integration involved in the time domain is quite cumbersome. Bogdan- 




O.l 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

PERCENTAGE ERROR, t 



Fig. (» — Convergence of mean-square response of processes ii(i) and x 2 (0- 
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off, and others, have shown that 

E[yl(t)} = ~ exp (-2W) £ ^# (20) 

V i x» 

where 7 ,-(2) and x> are both functions of a, , co, , and p. 

Sometimes it is more convenient to calculate the mean-square re- 
sponse by integration in the frequency domain than in the time domain. 
For example, Brown 17 has generalized Miller's 18 result for the process 
x s (t) to the case where $(£) is bounded on (— 00, 00) and integrable on 
every finite subinterval of (— 00 , 00). His frequency formulation is 

Cov v5 ft , L) = ± f Bfa, «,)£*(«, UiSA/c) do: (21) 



E[yl{t)} = A f \[ ff(«J*(«i - « 2 ) exp (i^t) d Ul 

07T J -~> J-m 



S,M du 2 (22) 



where H(io), $(w), and B(u, t) are the Fourier transforms of h(l), $(1), 
and h(t — t)4>(t), respectively, and the star denotes the complex con- 
jugate. Barnoski and Maurer, 19 using equation (21), evaluated E[y\(t)} 
numerically for cases where <f>(t) is the unit step and rectangular functions 
and f(t) is both white noise and noise with an exponentially decaying 
harmonic correlation function, that is, R { {t) = A exp ( — | r |) cos br. 
It was shown for white noise modulated by a unit step function that 
the system mean-square response will not exceed its stationary value 
for white noise. For correlated noise modulated in this same way, the 
system mean-square response may overshoot its stationary value. Note 
from equations (13) and (21) that 

0o (t, t') = Cov v5 (t, t'). 
Substitution of this relation into equation (12) yields 

E[yl(t)] = [ [ h(t - r)h(t - t') Cov„ 5 (t, r') dr dr' 

= -2 / /exp (fw«(T + t')) sin p(t — t) sin p(t — t') 

p Jr J r 

•Cov„ (r, t') dr dr*. (23) 

Substitution of g n in equation (13) into equation (12) yields 

E[y%t)} = A exp (-2ft*9 [ *(f) exp (2^) dr 
zp Jr 

- J I(t) exp (2fco n r) cos p(t - t) dr. (24) 
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Using the relation g a (r, t') = Cov„ 7 (r, /) and letting h (t) with param- 
eters w , to , and p be of the same form as h(t) in equation (11), it 
can be shown that 



E[yl(t)] 
1 



2p% 



exp (— 2fa„t) 



■j j | exp [($•«„ - focooXr + r')] sin p(« - t) sinp(* - r') 
• cospo(r - r') J 1(d) exp (2f o "o0) dd 

- j 1(6) exp (2r o co o 0) cosp (r + r' - 20) rf0 j dr dr'. (25) 

It should be noted that, in equations (24) and (25), I(t) = E[a 2 '(t)]j8(t) . 
The mean-square response to linear first-order moving average and 
autoregressive processes (assuming a sampling interval At) for first- 
order and second-order filters are also found. Let the transfer function 
for the first-order filter be h(t) = exp ( — a t) (corresponding to a dif- 
ferential operator p = d/dt + a ) and that for the second-order system 
be as given in equation (9), we obtain for the first-order ma process 



E[y 2 (t)] = ■£- [(1 + 2 ) - s*J(l - exp (-2a t)), (26) 



E[y\t)] = ^o exp (-2&J) 



(first-order filter); 



(1 + el)\~ (exp (2fco„0 - 1) - | exp (2&J) 

+ f cos 2p£ - (1 - f 2 )* sin 2p«J 

+ d l (1 - f 2 )* sin 2pt - f cos 2p< 

+ r exp (2f«„0 + i (1 - exp (2r«.*))J} , (27) 

(second-order filter); 
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and for first-order ar process 

fu,-7ai 1 [1 ~ (0! exp (ao))"'][(exp (-Oq)0i)' - 1] , 0fi . 

E \y (01 = f^ (In^Y-al ' (28) 

(first-order filter); 

7?n,,VAi - [exp (— o 1 Q(— a t sin pt - p cos a t t) + p] 

L u wj - (1 _ 0?)p . (a? + rtW + p>) 



[ exp (— a 2 t)(— a 2 si n pt — p cos a2<) + p] 
(1 - tf)pV + p 2 )(al + p 2 ) 

(second-order filter). 



(29) 



In the above 

c = AtE[g 2 ], 

q = 2 exp (a At), 

ai = fa>» + In & , 

^2 = f«n — In 0! . 

It follows from equations (26) and (27) that the mean-square response 
to the first-order ma process approaches a steady-state value 

lim E[y 2 (t)] = £- (1 - 0,) 2 , (first-order filter) ; (30) 

(-•oo 

= C ^ 4 ~t g,) , (second-order filter). (31) 

Similarly from equations (28) and (29) 

J& Em)] " d-^W-dn^ ■ " ° < * « 1. < 32 > 



I-»og 



(first-order filter); 

1 



(1 - tfMKln*,)' + <o 2 ] 2 - 4f 2 a, 2 (/n^) 2 } ' 
(second-order filter). 



(33) 



The above results indicate that the mean-square response of linear 
systems to either of the ar and ma processes approaches a certain sta- 
tionary value when the time of passage is sufficiently long. 
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IV. CONCLUSIONS 

When selecting a stochastic model it is important to consider not 
only the matching of statistical characterizations of real-world random 
data to those of the model, but also the effects of the model on vulner- 
able structural systems. The statistical characterizations are provided 
by joint probability distribution functions of the process conditioned on 
its duration and intensity. Ordinary power spectra and time-variable 
spectra such as the running spectrum 

S(ico, t) = I x(t) exp (—iur) dr 

Jo 

or the instantaneous spectrum 

p( W> = ^1 s(™, o r 

can be used to characterize and form the simulation basis for the sta- 
tionary and nonstationary processes, respectively. For the latter case, 
a second variable enters into the spectrum formulations and therefore 
complicates the analysis considerably. One alternative for analyzing 
nonstationary processes is to follow Priestley's theory of evolutionary 
processes and spectra. 20 " 22 Its applicability to ground motion and 
earthquakes will be investigated in a separate report. In general one 
must exercise engineering judgment upon consideration of the specific 
problem he studies in making an intelligent choice among all possible 
models. At the present time it appears that an earthquake ground- 
motion accelerogram can reasonably be treated as sectionally stationary 
when broken into three distinct phases. Each phase of the motion can 
be regarded as a short process and the corresponding power spectral 
density estimated by standard approaches to form the simulation bases 
for the x 2 {t) model. The applicability of this procedure is illustrated in 
Fig. 7 which shows the power spectrum densities associated with three 
distinct phases of recorded S21W ground motion during the Taft, 
California, 1952 earthquake. It is apparent from this earthquake that 
most of the input power to structures is provided by the midsection 
(from 3.3 to 13.6 seconds) of the motion. The power contained in low- 
intensity fluctuations preceding and following this stationary portion is 
relatively small. 

Finally it should be pointed out that the stochastic models investigated 
in this paper can be easily realized by using a computer. The response 
statistics of these models are also reasonably easy to find. Although the 
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Fig. 7 — Power spectral densities of Taft, California, S21W, July 21, 1952 earthquake. 

current study is concerned primarily with ground-motion simulation, 
the results obtained can be applied also to many other engineering 
problems when time series modeling is required. 
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